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The  eigenvalues  of  a  small  cluster  in  the  Cu03  plane  are  obtained  by  exactly  diagonalizing 
'-he  extended  Hubbard  Hamiltonian.  This  eigenenergy  spectrum  is  used  for  me  grand  parti¬ 
tion  function.  ~he  magnetic  and  thermodynamic  quantities  such  as  the  magnetic 
susceptibility  and  specific  heat  are  calculated  with  different  occupation  numbers  anu  dis¬ 
cussed  in  terms  of  the  excited  states  of  the  system. 
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I  .  Introduction 


Since  the  discovery  of  superconductivity  in  oxide  compounds  [1],  many  theoretical  and  ex¬ 
perimental  studies  were  carried  out  for  the  pairing  mechanism  on  these  compounds.  For 
most  crystalline  solids,  the  electronic  properties  are  well  described  by  the  band  theory  in 
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terms  of  a  one-body  Hamiltonian.  But  it  is  now  generally  believed  that  the  superconducting 
oxide  compound  do  not  lit  this  description.  In  this  case,  tne  excitations  of  these  systems  are 
totally  beyond  the  description  of  the  band-structure  approximation.  Therefore,  in  this  paper, 
we  study  the  magnetic  and  thermodynamic  properties  in  a  small  cluster  with  an  appropriate 
extended  Hubbard  Hamiltonian,  because  tnose  properties  mainly  depend  on  me  excitation 
spectrum  of  the  system. 

ttven  though  the  high-temperature  oxide  superconductors  have  some  similarities  with  con¬ 
ventional  superconductors,  there  seems  to  be  a  fundamental  difference  for  pairing  mecna- 
msm.  In  conventional  superconductors,  tne  attraction  between  two  electrons  arises  from  the 
phonon  exchange,  where  the  retarded  pnonon  attraction  is  able  to  overcome  the  screened 
Coulomb  repulsion  between  two  electrons.  However  it  is  very  hard  to  consider  a  similar 
mechanism  for  high-temperature  oxide  superconductors  with  short  coherent  length  c«  —  d. 
where  d  is  the  interparticie  distance  between  electrons.  The  wide  range  of  energy  gap  in¬ 
stead  of  a  sharp  gap  observed  in  most  tunneling  measurements  [2]  seems  due  to  the  fluctua¬ 
tion  effect,  which  is  expected  to  be  considerably  large  due  to  the  short  coherent  length.  More 
recently  discovered  n-type  superconductors[3],  which  have  similar  crystal  structures  to  the 
p-type,  provide  additional  information  about  these  compounds!  (1)  the  superconducting 
pairs  occur  in  purely  two  dimensions,  and  (2)  there  must  be  electron-hole  symmetry  for 
whatever  band  is  used  in  the  Hamiltonian. 

2  .  Theory 

i  he  central  physics  in  the  theory  oi  high-temperature  superconductivity  centers  on  which 
effective  Hamiltonian  one  should  take  to  describe  the  electrons  in  the  CuOj  plane.  The  high- 
Tc  compounds  can  be  explained  well  by  the  extended  Hubbard  Hamiltonian,  in  previous  pa- 
pers[4],  we  have  pointed  out  two  important  physical  features  of  the  pairing  of  electrons  and 
holes  in  n-type  and  p-type  oxide  superconductors,  respectively :( i )  me  vacuum  level  must 
„e  defined  as  (d’p4)  for  the  copper  ana  oxygen  sues  instead  of  the  closed  orbit  configuration 
id  p  )  used  Ly  other  l5],  x-ray  absorption  spectra  show  that  the  cnarge  carriers  i electrons ) 
are  mostly  located  on  the  copper  sues  <Cu  ' ;  in  n-type  superconductors[6],  while  the  charge 


earners  tholesi  are  tocatea  on  oxygen  su.es  tOl"i  in  p-iype  superconductors  i_7],  3y  this  ob¬ 
servation,  it  is  appropriate  to  define  me  vacuum  ievei  to  be  td9p‘),  wmch  is  an  unaopea  con¬ 
figuration  for  both  types  and  plays  a  critical  role  for  fixing  the  electron  ana  the  hole  energy 
levels.  (2)  Since  noles  are  not  true  particies  but  are  the  empty  states  which  represent  the 
aosence  of  electrons,  we  nave  to  be  more  careful  in  setting  up  the  Hubbara  Hamiltonian. 
Generally  the  on-site  Coulomb  repulsion  is  introduced  in  the  Hubbara  Hamiitoruan  wnen  two 
electrons  occupy  tne  same  orbit.  This  cannot  be  the  case  for  the  hole  states,  if  one  electron  is 
removed  from  two  electron  states  of  the  given  orbit,  i.  e..  one  hole  is  created,  the  on-site 
Eouiomb  repulsion  aisappears. 

.n  other  words,  me  on-site  CouiomD  repulsion,  wmch  is  usuaily  me  largest  parameter  in 
me  Hubbard  Hamiltonian,  appears  wnenever  a  hoie  is  created  from  the  given  closed  sneil 
orbit  instead  of  when  two  holes  are  created. 

Considering  the  above  two  features,  the  extended  Hubbard  Hamiltonian  can  be  written  as 

1  4  J  U  —  -  VU  t  H  It*  v.  p  X  n  , 
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where  t^  is  the  hopping  integral  between  neighboring  Cu  and  0  sites,  and  the  d.o's  ana  p,o 


are  second  quantization  operators  related  to  the  Cu  ana  0  sues,  respectively.  The  s', re¬ 
scripts  h  and  e  on  tne  Cu  site  operators  represent  the  holes  CuJ*(d’)  and  electrons  C  .  '  ( d:0), 
explicitly.  The  prime  on  a  summation  indicates  i  and  <i,  j>  indicates  nearen  neighbors 
between  the  Cu  and  0  sues.  Ed  and  E„  are  the  diagonal  kinetic  energies  of  the  id  anu  bp  or¬ 
bitals,  respecitvelv.  The  eiectron-hole  pair  creation  potential  L'*  in  the  last  term  of  the  Ham¬ 
iltonian  is  introduced  to  account  for  the  largest  valence  chanre  dif f  .once  wnen  an  electron 


>n  an  0  site  moves  to  a  neighboring  Cu  sue.  ie.  e„  Cu;'fd’)  ~;p’)—  Cu'*M'0'  -■  '  pH. 

As  mentioned  earner,  the  on-site  Coulomo  repulsions  L'd  ar  .  L'p  appear  when  an  electron  is 
created  on  a  Cu  sue  or  a  hole  is  created  on  an  0  sue.  Wi.h  the  condition  E„-  l'.,  >E,-?  V  the 
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undoped  LajCuO<  and  NbjCuO*  nominally  consist  of  Cu:"and  0J~.  and  the  holes  introduced 
by  doping  with  the  divalent  Sr  in  La2CuO«  occupy  the  oxygen  orbitals  in  p-type  superconduc¬ 
tors,  while  the  electrons  introduced  by  doping  with  the  tetravalent  Th  in  NdjCuO*  occupy  the 
copper  orbitals  in  n-type  superconductors.  Clearly  this  observation  is  due  to  the  presence  or 
absence  of  electrons  in  the  antibonding  states  of  Cu(d^-r2)  and  Ofp,  or  pT)  orbitals. 

We  believe  this  hybridization  is  strongly  related  to  the  quick  disappearance  of  the 
antiferromagnetic  background  on  the  Cu  sites  with  a  small  amount  of  doping.  Even  though 
the  direct  hopping  integral  (t«) between  neighboring  oxygen  sites  is  not  negligible,  the  direct 
oxygen  band  is  supposed  to  be  much  lower  than  the  antibonding  states  of  the  Cu  and  0  orbit¬ 
als^]  due  to  the  condition  E<-r  L’4>E,-r  LV  Therefore,  the  direct  oxygen  hopping  integral  is 
not  included  in  the  Hamiltonian. 
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Fig  I.The  cluster  considered  in  the  calculation  includes  two  unit  cells  in  the  Cu02  square  lat¬ 
tice.  The  copper  sites  and  oxygen  sites  are  represented  by  the  circles  and  plus  signs, 
respectively. 

Based  on  the  above  arguments,  we  have  proposed  a  local  pairing  mechanism  and  estimat¬ 
ed  the  binding  energy  of  two  carriers  in  a  previous  paper[4].  In  the  present  study,  we  exam¬ 
ine  the  magnetic  and  thermodynamic  quantities  such  as  the  magnetic  susceptibility  and  spe¬ 
cific  heat  by  exact  diagonalization  on  small  two-dimensional  clusters  of  two  unit  ceils  on  a 
CuOj  lattice  as  shown  in  Fig.  1.  By  introducing  the  chemical  potential  u  into  the  Hamiltom- 
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an,  we  can  ailow  the  numoer  ol  electrons  to  vary  in  this  system.  In  this  case,  ah  the  second 
quantization  operators  in  Eq.(l)  can  be  changed  to  the  true  electron  representation  without 
the  hole  picture,  information  on  the  interaction  between  the  cluster  and  the  rest  of  the 
system  is  also  included  in  the  chemical  potential.  The  Hubbard  Hamiltonian.  Eq.(l),  of  this 
system  is  now  modified  in  the  presence  of  a  uniform  magnetic  field  as 
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where  r  is  the  magnetic  moment  of  the  electron,  and  the  uniform  magnetic  field  B  is  perpen¬ 
dicular  to  the  CuO,  plane. 

The  eigenvalues  of  the  above  Hamiltonian  can  be  obtained  by  diagonalizing  exactly  in  a 
small  cluster  system.  The  grand  partition  function  Z  can  be  written  immediately  with  the 
eigenvalues.  To  calculate  tne  magnetic  and  thermodynamic  quantities,  we  have  to  first  evalu¬ 
ate  the  chemical  potential  of  the  system  tor  the  given  occupation  numbers.  The  chemical  po¬ 
tential  can  be  determinea  by  inverting  the  expression 


N)  = 


K0T 


r  d  tnZ 
-  dp 


T.  V 


(3) 


in  terms  of  the  occupation  number  N  and  the  other  parameters  in  the  Hamiltonian.  The  mag¬ 
netic  susceptibility  x  and  the  specific  heat  C.  can  be  calculated  with  a  knowledge  of  the 
grand  partition  function  and  the  chemical  potential  in  the  usual  manner. 


ffl .  Numerical  Calculations  and  Discussion 

Considering  four  possible  states  per  site,  the  given  clusters  has  4N  eigenvalues  and 
eigenvectors.  Since  the  Hamiltonian  is  isotropic  in  spin  space,  the  total  spin  S  and  its  z-com- 
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ponem  Sz  are  good  auantum  numbers.  Then  the  4s  x  4*  Hamiltonian  matrix  can  De  reaucea 
to  many  blocks  with  different  number  of  occupations.  In  the  following  calculations,  we  nave 
used  these  reported  parameters,  which  are  in  the  most  acceptable  range  [8]:  aE  =  1.5  eV.  •, 


„  =  1.2  eV,  Ud=  9-0  eV.  L'„  =  6.0  eV  and  two  different  values  of  (1)  U«=  O.OeV.  (2)  1'*  = 
1.5  eV.  In  all  cases,  it  is  found  that  the  ground  state  in  the  state  is  the  state  with  minimal 
Sz  |  .  This  result  is  due  to  the  relatively  large  hopping  integral  t*  and  more  degree  of  free¬ 
dom  with  less  S,  I  .  The  behavior  of  the  magnetic  and  thermodynamic  quantities  in  the 
cluster  model  can  be  understood  from  the  excitation  spectra.  In  Figs.  2  and  3.  we  have 
shown  the  magnetic  susceptibility  and  specific  heat  curves  of  the  two  unu-ceil  model  for  N  = 
10,  which  is  a  vacuum  cofiguration  in  our  aefiniuon.  The  low-temperature  peaxs  in  these  fig¬ 
ures  are  due  to  the  transition  from  tne  ground  state  of  singlet  spin  configuration  on  two  cop¬ 
per  sites  to  the  triplet  spin  configuration  wmcn  is  the  first  excited  state.  By  contrast,  tnese 
peaks  disappear  for  N=  9  and  N=  6,  which  are  shown  in  Figs.  4  through  7. 
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Rg  2.  Magnate  auaceptibCty  arves  (r/2Fl  for  N=  10  with  (1)  L',*  =  O.OeV  and  (2)  12*  = 
1.5eV. 


-84- 


-00  COG  SCO  •  :co 

[MPERATURE  (k) 


Pg.  3  Specific  heat  arve  (Cv/N)  for  N=  10  with  ( 1 )  U,*=  O.OeV  and  (2)  U*  =  1.5eV. 
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:,Q-  4  Magnetic  susceptJbiStv  cxrves  ( x/2r)  for  N=  9  wuh  ( 1)  t\»  =  0.0  eV  and 
1.5eV. 
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Fig.  7.  Specific  heat  curv«*  (Cv/N)  for  N=  8  with  (1)  U*=  O.OeV  and  (2)  U*=  1.5eV. 


For  more  than  half-filling  cases  (N>6),  the  energy  spectrum  can  roughly  be  divided  into 
some  groups  with  spacings  of  on-site  Coulomb  repulsion  Ud  or  Up,  and  the  spacings  between 
successive  levels  within  the  same  group  are  expected  to  be  on  the  order  of  Upd/Ud  or  i‘pd/  Up. 
As  it  turns  out  in  our  numerical  calculations,  the  energy  spacings  between  the  ground  state 
and  the  first  excited  state  are  greater  than  0.5  eV  for  all  cases  except  the  vacunm  comigura- 
Mon  N=  10.  The  energy  difference  between  the  ground  state  (S,=  0)  and  first-excited  state 
S,=  i)  results  in  30  meV  — 60meV  for  N=  10  with  various  parameters  in  the  Hamiitoman. 

This  energy  difference  gives  the  superexchange  parameter  J  between  neighboring  copper 
sites  in  the  Heisenberg  representation, 

H=  J£  Si, -5/  (4) 

(c  J) 

which  is  a  good  approximation  for  an  undoped  matenai[9].  J  can  be  obtained  alternatively 
by  the  fourth-order  perturbation  approximation  to  Eq.(l)[4], 
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which  has  very  good  agreement  wnn  me  aoove  numerical  result. 

in  Fig.  4.  two  magnetic  suscepiiDiiiiy  curves  almost  overlap  ana  cecreases  monotomcaiiy 
as  T  .ncreases  because  the  grouna  ievei  nas  i  minimal  Sz=  1  2.  which  is  aifferent  :rom  me 
case  of  the  N  =  8  and  N  =  10.  In  Fig.  6.  tne  magnetic  susceptibility  are  increasing  wu.i  in¬ 
creasing  temperature  as  in  Fig.  2.  but  nc  peaxs  appear  unaer  1000  k  aue  to  the  large  excita¬ 
tion  energy  from  the  ground  state,  as  explained  earlier.  These  qualitatively  different  types  of 
behavior  of  the  magnetic  susceptibility  are  merely  caused  by  the  smaiiness  ot  me  <  .  ter. 
With  the  same  explanation  of  large  excitaiton  energy  from  the  grouna  state,  tne  specific 
neat  curves  in  Figs.  5  and  7  have  no  Deatcs  in  me  iow-temperature  regime. 

IV .  Conclusions 

The  extended  HuDbard  Hamiltonian  nas  oeen  solved  for  ciusters  of  two  unit  cells  by  an 
exact  diagonalizauon.  The  ground  states  are  always  related  to  the  minimal  i  Sz  i  configura¬ 
tion.  The  energy  spectrum  is  grouped  with  the  order  of  on-site  coulomb  repulsion  U„  or  Ua. 
and  the  successive  levels  in  each  group  are  separated  at  the  order  of  tWCd  or  t^'LV  But 
the  first  excited  state  in  an  undoped  regime  is  the  order  of  superexchange  energy  between 
nearest-neighbor  copper  sites,  which  agrees  well  with  the  result  obtained  by  the  fourth-order 
perturbation  calculation.  This  small  excitation  energy  contributes  to  the  low-tempertaure 
peaxs  in  tne  magnetic  suscepuDiiity  ana  scecific  neat  curves  tor  N  =  10.  Due  to  tne  .arge  ex¬ 
citation  energy  tor  tne  other  configuration,  no  peaxs  appear  in  the  magnetic  susceptibility 
and  specific  heat  curves  in  the  low-temperature  regime.  From  these  resuits,  we  may  con¬ 
clude  that  the  iow-temperature  phyr.  ai  anomalies  appear  oniy  near  the  unaopea  regime. 
Thu  magnetic  susceptibility  curve  for  the  case  of  N  =  9  is  essentially  of  the  Curie  form, 
monotomcaiiy  decreasing  as  tempenure  increases.  This  is  due  to  the  odd  number  oi  electrons 
in  a  small  cluster.  We  believe  that  this  kind  of  behavior  will  be  diminshed  if  we  increase  the 
size  of  the  cluster. 

Finally,  we  would  sav  that  the  results  ootained  in  this  calculation  are  all  very'  stable  with  a 
wide  range  of  parameters  in  the  Hamiltonian. 
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Therefore,  the  qualitative  nature  of  the  results  does  not  change  within  tne  acceptaole  pa¬ 
rameter  ranges,  even  though  we  have  snown  oniy  two  different  values  of  L,,  wnich  is  sup¬ 
posed  to  be  the  most  uncertain. 
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